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I. INTRODUCTION 



Life at low Reynolds numbers has attracted considerable attention in the past few 
decades 1-4 . The absence of inertia plays a remarkable role in the swimming of microor- 
ganisms. In a Newtonian flow, the scallop theorem 1 constrains the types of locomotion 
strategies which are effective in the microscopic world and reciprocal motion - as are called 
those with a time-reversal symmetry - cannot lead to any net propulsion (or fluid transport). 
Microorganisms evolved different propulsion strategies to achieve micro-propulsion, includ- 
ing the active propagation of flagellar waves for eukaryotic cells and the passive rotation of 
rigid helical flagella for bacteria 2 . 

The physics of low Reynolds number locomotion is relatively well explored in the New- 
tonian limit (see reviews in Refs. 1-4 and references therein). Beyond improving our under- 
standing of biological processes, applications of these physical principles led to progress in 
the design of synthetic micro- swimmers for potential future biomedical applications 4-8 . In 
contrast, fundamental properties of life in complex, non- Newtonian, flows remain surpris- 
ingly unaddressed. Non-Newtonian flow behaviors can be appreciated through well-known 
manifestations from daily life, for example the climbing of dough up kitchen mixing blades 
(termed rod-climbing, or Weissenberg, effect) or the remarkable behavior of Silly Putty, a 
popular toy which bounces like a solid rubber ball when thrown to the floor but melts like 
a fluid when left on a surface for some time 9-11 . 

Many situations exist wherein microorganisms encounter biological fluids which have 
polymeric microstructures and non-Newtonian rheological properties. For example, sperma- 
tozoa swim through the viscoelastic cervical mucus and along the mucus-covered fallopian 
tubes 3 ' 12-16 ; cilia lie in a layer of mucus along the human respiratory tract 17 ; Helicobacter 
pylori, a bacterium causing ulcer, locomotes through mucus lining of the stomach 18 ; spiro- 
chetes moves through host tissue during infection 19 ; in biofilms, bacteria are embedded in 
cross-linked polymer gels 20-24 . 

Physically and mathematically, the presence of polymeric stresses in a complex fluid 
means that the usual properties associated with the absence of inertia in the Newtonian limit 
cease to be valid, in particular kinematic reversibility and the linearity of the flow equations. 
In return, non-Newtonian effects such as stress relaxation, normal stress differences, and 
shear-rate dependent viscosity manifest themselves 9-11,25 . 
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Past theoretical and experimental studies have investigated the waveforms and swimming 
paths of microorganisms in complex fluids 3,4 ' 12-15 ' 26-28 . An active discussion in the biome- 
chanics community has recently focused on the simple question: does fluid elasticity enhance 
or deteriorate propulsion at the microscopic scale? Theoretical studies on infinite models 29-31 
showed that, for fixed body-frame kinematics, the propulsion speed decreases in a viscoelas- 
tic fluid. Numerical studies on a finite swimmer 32 demonstrated that the propulsion speed 
could be enhanced by the presence of polymeric stress for some prescribed kinematics. Ex- 
perimental investigations suggested evidences for both 33,34 . It was also shown that reciprocal 
actuation on a fluid, unable to provide net locomotion or flow transport in the Newtonian 
case, can be rendered effective by viscoelasticity 35-37 . The presence of polymeric stress has 
also interesting consequences on the rate of flagellar synchronization 38 . 

Normal stress differences in a complex fluid are responsible for a number of important 
non-Newtonian effects 9-11 including the rob-climbing effect mentioned above impacting many 
applications such as mixing, and the swelling of polymer melts when extruded from dies 
in manufacturing processes posing constraints on the rate of extrusion. In a pure shear 
flow with an arbitrary Reynolds number and a shear rate 7, assuming that the flow is in 
the x— direction and the velocity varies in the y— direction, the z— direction being called 
the neutral direction 9 , the first and second normal stress coefficients are defined as ^1 = 
{jxx — T yy)/i 2 an d ^2 = {jyy — T~ zz )/i 2 respectively, where r^- are the components of the 
deviatoric stress tensor. In a Newtonian flow, there are no normal stress differences (^1 = 
\&2 = 0), whereas for polymeric fluids typically ^ > and ^ 2 < 0. The magnitude of 
the second normal stress coefficient is usually much smaller than that of the first normal 
stress coefficient (l^l ^i). In the rob-climbing phenomenon, both first and second 
normal stress coefficients contribute to the effect 9 . However, due to its small magnitude, the 
effect of the second normal stress coefficient is shadowed by that of the first normal stress 
difference 9 . The existence of the second normal stress difference can be demonstrated in 
a free-surface flow driven by gravity through a tilted trough: a Newtonian fluid has a flat 
free surface (with negligible meniscus effect), while the free surface of a non-Newtonian fluid 
becomes convex due to second normal stresses 39-41 . 

In this work, we propose a simple mechanism able to take advantage of the presence 
of normal stress differences to propel in a complex fluid. Our geometry, shown in Fig. 1, 
consists of two linked small spheres propelling under the action of an external torque, a setup 
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we will refer to as a "snowman". Locomotion is enabled solely by the presence of normal 
stress differences, and no motion exists in a Newtonian environment, a fact that can in turn 
be used to infer the normal stress coefficients of a complex fluid. In essence, as complex 
fluids lead to new modes of small-scale propulsion, symmetrically the presence of propulsion 
in an environment can be used to locally probe the rheological properties of the fluid. 

The first normal stress coefficient of a fluid can be measured directly from a conventional 
cone-and-plate rheo meter 9-11 ; the measurement of the second normal stress coefficient how- 
ever has been a longstanding challenge 42-45 . A number of methods were proposed (see a re- 
view in Ref. 45 ), including a modified cone-and-plate rheometry with pressure transducers 9 ' 43 , 
a subtle evaluation of a combination of cone-and-plate and parallel-plate experiments 42 , 
rheo-optical measurements 42 ' 44 , and the use of a cone-and-partitioned plate tool 45 . Re- 
cently, a microrheological technique was proposed to measure the first and second normal 
stress coefficients 46 . In microrheology, colloidal probes are either actively driven, or passively 
diffusing, and their dynamics allows to infer local rheological information. Microrheology 
enjoys many advantages over conventional macroscopic rheological measurements 47 ' 48 , in- 
cluding the reduction in sample size, the ability to probe spatially-inhomogeneous environ- 
ments, and the possibility of performing measurements in living cells 47-50 . The mechanism 
we propose in this paper would be classified as "active" microrheology, a situation where 
colloidal probes are actively manipulated to drive the material out of equilibrium and probe 
its nonlinear mechanical properties 47 ' 51 . We offer in this paper an alternative microrheo- 
logical technique capable of probing both first and second stress coefficients by using only 
kinematic measurements. 

This paper is organized as follows. In Sec. II, we introduce the geometric and kinematic 
setup of our proposed mechanism and the polymeric fluid models adopted in our study. 
We first investigate in Sec. Ill the propulsion characteristics of the snowman in a complex 
fluid, followed in Sec. IV by the method of inferring the normal stress coefficients from its 
locomotion. We then provide a qualitative, and intuitive, explanation of the locomotion 
enabled by normal stresses in Sec. V before concluding the paper in Sec. VI. 
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FIG. 1. Geometrical setup of two spheres ("snowman") rotating with angular velocity f2 along their 
separation axis. The radii of the upper and lower spheres are denoted by Rjj and Rl respectively. 
The centers of the spheres are separated by a distance, h (for touching spheres, h = Ru + Rl)- 



II. SETUP 
A. Kinematics 



By symmetry, the rotation of a single sphere in any homogeneous fluid produces no net 
locomotion. Inserting a second sphere, of different size, breaks the geometrical symmetry 
and can potentially allow locomotion. We first consider the rotation of two unequal spheres 
touching each other as a single rigid body (see the geometry and notations in Fig. 1), the 
"snowman" geometry. We label the line of centers of the spheres as the z-axis. Without 
loss of generality, we assume the radius of the upper sphere (Ru) is smaller than that of the 
lower sphere (Rl > Ru)- The distance between the centers of the spheres is denoted by h. 
For the case of touching spheres, we thus have h = Ru + Rl- 

From a kinematic standpoint we assume that the rigid body rotates with a steady angular 
velocity about the z-axis, Q = (0, 0, Q > 0), but is otherwise free to move. Given that the 
snowman is axisymmetric, the only direction it could potentially move is the ^-direction. 
We assume the net hydrodynamic force acting on the snowman is zero for all times (free- 
swimming condition), and aim at computing the rigid body (swimming) velocity necessary 
to maintain force-free motion. 

In a Newtonian fluid without inertia, it is straightforward to show using kinematic re- 
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versibility and reflection symmetry that a rotating snowman cannot swim - a result true for 
any degree of geometrical asymmetry. The central question at the heart of this paper is: 
Can elasticity of the fluid enable propulsion of the snowman? We answer this question in 
the following sections by studying the locomotion of a rotating snowman in polymeric fluids 
described by the two constitutive relations. 

B. Polymeric fluid dynamics 

We consider an incompressible low- Reynolds-number flow in a complex fluid. Denoting 
the velocity field as u and the fluid stress as a = —pi + r, where p is the pressure, and 
t is the deviatoric stress tensor, the conservation of mass and momentum are given by the 
continuity equation and Cauchy's equation of motion respectively 

V-u = 0, (1) 
V • cr = 0. (2) 

For closure, we require a constitutive equation relating the deviatoric stresses r to the 
kinematics of the flow. Obviously a large number of models have been proposed in the past 
to describe polymeric fluids. In this work two constitutive equations are used to study the 
viscoelastic locomotion of a snowman. 

1. Oldroyd-B fluid 

The classical Oldroyd-B constitutive equation is arguably the most famous constitutive 
model for polymeric fluids 9-11 ' 25 . It has a sound physical origin and can be derived from 
a kinetic theory of polymers in the dilute limit by modeling polymeric molecules as lin- 
early elastic dumbbells. The predictions also agree well with experimental measurements 
up to order one Weissenberg numbers, although it is known to suffer deficiencies for larger 
values 9-11,25 . In an Oldroyd-B fluid, the deviatoric stress is the sum of two components, 
t = t s + t p , where r s and r p denote, respectively, the Newtonian solvent contribution and 
polymeric contribution to the stress. The constitutive relation for the Newtonian contribu- 
tion is given by r s = rjs'j, where 7 = Vu + Vu T is the rate of strain tensor and T] s is the 
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solvent contribution to the viscosity. The momentum equation can thus be written as 

-Vp + ri s V -7 + V- r p = 0, (3) 

The polymeric stress t p is then assumed to be governed by the upper-convected Maxwell 
equation 

t p + A r p = Vp <y, (4) 

where A is the polymeric relaxation time and t] p is the polymer contribution to the 

viscosity 9-11 . The upper-convected derivative for a tensor A is defined as 

v 3 A 

A= ^ + u- VA-(Vu T - A + A- Vu), (5) 

ot 

which calculates the rate of change of A while translating and deforming with the fluid. 

Combining the Newtonian and polymeric constitutive relations, we obtain the Oldroyd-B 
constitutive equation for the total stress, r, as 

v v 

T + \ T= 7](j + \ 2 j), (6) 

where the total viscosity is given by rj = rj s + t] p , and A2 = \( denote the retardation times 
(we define the relative viscosity ( = r] s /r] < 1). For steady shear of an Oldroyd-B fluid, 
both the viscosity and the first normal stress coefficient are constant, and the second normal 
stress coefficient is zero 9 . The Oldroyd-B fluid is the model we will use for our numerical 
approach. 



2. Second- order fluid 

For slow and slowly varying flows, the second-order fluid model applies. It is the first 
non-Newtonian term in a systematic asymptotic expansion of the relationship between the 
stress and the rate of strain tensors called the retarded-motion expansion. It describes small 
departures from Newtonian fluid behavior, and the instantaneous constitutive equation is 
given in this model by 

1 v 

r = 777- ~# l7 + #3(7.7), (7) 

where #1 and #2 are the first and second normal stress coefficients respectively. Note that 
if A = while A2 7^ in the Oldroyd-B model, Eq. (6), it reduces to a second-order fluid 
with a vanishing second normal stress coefficient (#2 — 0) 9 . The second-order fluid model 
will enable us to derive theoretically the behavior of the snowman for small deformations. 
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C. Non-dimensionalization 



We non-dimensionalize lengths by the radius of the lower sphere Rl, times by and 
use the total fluid viscosity, r/, to provide the third fundamental unit. Hence, velocities, shear 
rates, and stresses are scaled by Rl&, and r)Q respectively. The dimensionless radius of 
the upper sphere becomes then r* = Rjj/Rl while the lower sphere has now radius 1. We 
have h* = h/R^ denoting the dimensionless distance between the centers of the sphere 
(h* = 1 + r* for two touching spheres). Both spheres rotate at the same dimensionless 
unit speed, Q* — 1. The starred variables represent dimensionless variables in this paper. 
The Deborah number 9-11 , De = XQ, is a dimensionless number defined as the ratio of a 
characteristic time scale of the fluid (the polymeric relaxation time, A) to a characteristic 
time scale of the flow system and appears in the dimensionless momentum equation 

and upper-convected Maxwell equation 

-Vp* + CV ■ 7* + V ■ t p * = 0, (8) 

r p * + Der Vp *= (1 -C)7*. (9) 

The limit De = corresponds to a Newtonian fluid. 

Alternatively, the upper-convected Maxwell equation of the polymeric stress, Eq. (9), 
can be combined with the constitutive relation of the Newtonian contribution to obtain the 
dimensionless Oldroyd-B constitutive equation for the total stress r* as 

r* + De r*= 7* + De 2 7*, (10) 

where we have defined another Deborah number, De2, in terms of the retardation time, 
De 2 = A 2 fi = De(. 

The dimensionless constitutive relation for a second-order fluid is now given by 

t* = Y- De so U* +Bj* ■ 7*) , (11) 

where we have defined another Deborah number for the second-order fluid, namely De so = 
and B = -2* 2 /*i > 0. 
Importantly, we note that the definition of the Deborah number of an Oldroyd-B fluid is 
different from that of a second-order fluid, because the relaxation time of an Oldroyd-B fluid 
is defined only by the polymer, whereas the relaxation time of a second-order is defined by 
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both the polymer and the solvent 52 . The two Deborah numbers are related by the relation 
De so = De(l — £)■ We shall mostly use the Deborah number defined for an Oldroyd-B fluid 
(De) for the presentation of our final results, since we feel it is the one with the most intuitive 
definition. The Oldroyd-B equation is valid up to moderate De, and the second order fluid 
is valid for small De (or De so ), and we thus expect the results from both models to match 
when De (or De so ) is sufficiently small. 

III. PROPULSION OF SNOWMAN IN A COMPLEX FLUID 

As argued in Sec. II A, asymmetry alone does not lead to net locomotion upon rotating 
a snowman in a Newtonian fluid. We now explore the effects of fluid elasticity on the 
propulsion of a snowman: Does it even move? Which direction does it go? And how fast? 
Using the Oldroyd-B fluid model, we first explore numerically the propulsion characteristics 
of the snowman from small to moderate Deborah numbers. Next, the small De limit is 
studied analytically via the second-order fluid model. 

A. Moderate Deborah number 

We employed a finite element model to compute the polymeric flow as described by 
Eqs. (8) and (9). A formulation called the Discrete Elastic- Viscous Split Stress (DEVSS- 
G) 53,54 is implemented here to improve numerical stability. The momentum equation, 
Eq. (8), is rewritten as 

V • /i a (Vu* + Vu* r ) - Vp* + V ■ - V • (/i a - C)(G + G T ) = 0, (12) 

where the tensor G = Vu* is introduced as a finite element approximation of the velocity 
gradient tensor Vu*. An additional elliptic term, V • fi a (Vu* + Vu* r ) - V • /u (G + G T ), 
is added into the momentum equation for stabilization 55 . In the limit that the mesh size 
in the finite element approximation tends to zero, G approaches Vu* and the elliptic term 
vanishes, reducing Eq. (12) to Eq. (8). G is also used to approximate the velocity gradient 
term Vu* in the constitutive equation, Eq. (9). For simulations in this work, we choose 
li a — 1 as in Liu et a/. 54 . 

A Galerkin method is used to discretize the momentum equations, continuity equation, 
and the equation for the additional unknown G. Quadratic elements are used for u* and lin- 
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ear elements for both p* and G. The streamline-upwind/Petrov-Galerkin(SUPG) 56 method 
is adopted to discretize the constitutive equation, Eq. (9), to improve numerical stability. 
The resulting weak form of the model is formulated as 

{S + -^u* ■ VS, r p * + De(u* • Vt p * - G T ■ t p * - t p * ■ G) - (1 - C )(G + G T )} = 0, 

(13) 

where S denotes the test function for t p *, h c is a characteristic mesh size, and U c is the 
magnitude of a local characteristic velocity (we choose the norm of u* as U c ). The frame- 
work for the implementation is provided by the commercial software COMSOL, which was 
successfully used for simulating the locomotion of squirmers in a viscoelastic fluid at low 
Reynolds numbers 57 . 

We perform three-dimensional axisymmetric simulations on a two-dimensional mesh con- 
structed with triangle elements. Sufficiently refined mesh is generated near rotating objects 
to resolve the thin stress boundary layers, necessary to overcome numerical instabilities 58 ' 59 
and improve accuracy. We validate our implementation by comparing numerical and ana- 
lytical values of the hydrodynamic torque on a rotating sphere in the Newtonian fluid. For 
the viscoelastic model, we validate our approach against the simulations in Lunsmann et 
al. m of a sedimenting sphere in a tube filled with Oldroyd-B fluid and the analytical results 
in Bird et al. 9 of a rotating sphere in a second-order fluid. 

Equipped with our computational model, we are able to show that fluid elasticity does 
indeed enable the propulsion of the snowman provided the two spheres have unequal sizes 
(r* < 1). The snowman always swim in the positive z-direction (see Fig. 1), i.e. from 
the larger to the smaller sphere. For illustration, we compute the dimensionless propulsion 
speed, U* = U/RlQ, of a typical snowman (r* = Rjj/R l = 0.5) as a function of the 
Deborah number, De (dot-dashed line - red online, Fig. 2), for a fixed relative viscosity 
( = 0.5. When De = 0, the fluid reduces to the Newtonian limit and we recover that no 
propulsion is possible in this case. For small values of De, the propulsion speed appears 
to grow linearly with De, a result confirmed analytically in the next section. A maximum 
swimming speed is reached at De ~ 1.75, before decaying as De continues to increase. 

In addition to the primary flow (the Newtonian component, De = 0), elastic stresses 
around the snowman generate a secondary flow, understood simply as the difference between 
the total flow and the Newtonian component. A typical secondary flow pattern is shown 
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FIG. 2. Demonstration of snowman locomotion. In the case r* = Ru/Ri = 0.5 and £ = 0.5, 
we plot: (a) Dimensionless propulsion speed, U/Rl$1, as a function of the Deborah number, De. 
Dot-dashed line (red online): numerical simulations in an Oldroyd-B fluid; solid line (blue online): 
theoretical calculation using the reciprocal theorem in a second-order fluid, Eq. 36; (b) The stream- 
line pattern and speed (shaded/color map) of the secondary flow for De = 0.1 (streamline patterns 
at higher De are qualitatively similar). 

in the frame of the snowman in Fig. 2b (De = 0.1 and ( = 0.5). We depict the velocity 
vectors and streamlines with the shaded/color map representing the flow speed. Fluid is 
drawn towards the snowman parallel to the equatorial plane and then expelled along the 
axis, while a ring vortex is detected in the front. The maximum speed of the secondary flow 
is observed at the rear of the snowman, only about 0.7% of the characteristic speed of the 
primary flow R^VL. 

B. Small Deborah number 

To provide a theoretical approach to the snowman locomotion and to quantify the con- 
nection between locomotion and rheology in the following sections we now consider the 
second-order fluid, which we remind is valid in the small-De limit only (Eq. 11). All vari- 
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ables are expanded in powers of the Deborah number, De so , as 



<j = cr + De so <Ti + 0(De so 2 ), 
u = u + De so ui + 0(De so 2 ), 
7 = 7 + Be so ^f 1 + 0(De so 2 ), 
U = U + De so U 1 + O(De so 2 ), 



(14) 
(15) 
(16) 
(17) 



where U denotes the propulsion velocity, U = (0,0, U). Other variables are expanded 
similarly. We drop the stars hereafter for simplicity, and all variables in this section are 
dimensionless unless otherwise stated. The locomotion problem is then solved order by 
order. 

1. Zeroth-order solution 

The zeroth order solution, {cr = — PqY + 7 , u }, satisfies the Stokes equations, 



where er = — p I + 7 . This is the Newtonian flow for two touching spheres rotating (at a 
rate of fi) about the line of their centers (z-axis). The exact solution in terms of analytical 
functions was given by Takagi 61 in tangent-sphere coordinates. No propulsion occurs in the 
Newtonian limit, U = 0, as expected. 

2. First-order solution 

The first order solution (cr^ ui) of the main problem satisfies 



V ■ <r = 0, 



V ■ u = 0, 



(18) 
(19) 



V • ui 











(20) 
(21) 



where 



v 



o-i = -pil + 7i - 7o - B io ■ 7o- 



(22) 
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To compute the value of the first order propulsion velocity, Ui, we will use a version of the 
reciprocal theorem for Stokes flows adapted to self-propulsion in viscoelastic fluids 36 ' 46 ' 62-68 . 
Consider an auxiliary problem with identical geometry, {cr aux , u aux }, satisfying 

V ■ <x aux = 0, (23) 

V ■ u aux = 0. (24) 

Taking the inner product of Eq. (20) with u aux , minus the inner product of Eq. (23) with 
ui, and integrating over the entire fluid volume, we have trivially 

/ u aux • (V • <ri) - ui • (V ■ <r^)dV = 0. (25) 
Jv f 

Using vector calculus we can rewrite the integral in the following form 69 

/ V ■ (u aux ■ <Ti - ui • (T aux )dV = / (Vu aux : (Ti - Vui : (T aux )dV. (26) 
JVf Jv f 

The left-hand side of Eq. (26) can be converted to a sum of surface integrals by the diver- 
gence theorem while the right-hand side can be simplified using the first-order constitutive 
equation, Eq. (22), leading to 



^2 J n • ( u aux • <Ti - Ui ■ Vau X )dS = J ^7 + .070 ■ 7 ^ : Vu 



dV, (27) 



where S a denotes the surface of different spheres (a = 1,2) and n represents the outward 
normal vector on the surface. The important simplification which took place in the right 
hand-side of Eq. (26) is that all Newtonian terms included in cr aux and <T\ have canceled 
each other out by symmetry, and thus the only piece remaining in the right-hand side of 
Eq. (27) is the non-Newtonian contribution 63 ' 68 . 

Now, let Ui and fii be the (unknown) first order translational and rotational velocities 
of the spheres in our main problem, while the translational and rotational velocities of the 
spheres in the auxiliary problem (known) are given by U aux and f2 aux . On the surface S a of 
one sphere, the no-slip and no-penetration boundary conditions lead to 

u aU x = U aux + f2 aux x r, (28) 
ui = Ui + fli x r, (29) 
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where r is the position vector describing the surface. The integral relation, Eq. (27), becomes 

[ n-ndS + O^- [ Tx(n-tn)dS 

q, J Set J S a 

-U? • / n ■ v aux dS - n? ■ / r x (n • <T aux )d,S 



7o + 5 7o • 7o ) :Vu 



aux 



tiV. (30) 

In Eq. (30), the integrals J s n-cridS and J 5 rx (n-<Ti)dS represent the net hydrodynamic 
force and torque acting on the sphere a by the first order flow field. Let us denote F" = 
— j s n • a\dS and T" = — j s r x (n ■ <Ji)dS the net external force (F") and external torque 
(T°) acting on each sphere; the appearance of a minus sign comes from the fact that the 
total force and torque (external + fluid) acting on a body have to sum to zero in the absence 
of inertia. In the free-swimming case there is an additional stronger constraint, namely the 
total external force (or equivalently, the total fluid force) has to remain zero at all instant 
(we will enforce this constraint shortly). Defining also F" ux and T" ux as the external force 
and torque required to balance the fluid drag and torque on each sphere in the auxiliary 
problem we see that Eq. (30) is transformed into 

-u^ ux ■ f? - • t? + u? ■ f^ ux + n? • t^ ux = jf ( 7o + b% • <y ) = Vu aux dv. 

(31) 

The above relation remains actually true for any number of spheres and kinematics. In 
the case of a snowman, we have two spheres connected as a rigid body in both the main 
and auxiliary problems, hence XJ\ = U 2 = Ui, CI 1 = Q 2 = f2i, U^ ux = U 2 UX = U aux , and 
f2 aux = f2 aux = ^aux- I n the main problem we impose a rotational rate fl on the snowman, 
which has been accounted for in the zeroth order (Newtonian) solution, hence ft" = for 
all n > 1. In addition, we define in the main problem the total external force and torque 
acting on the rigid body as Fi = F 1 + F 2 , Ti = T 1 + T 2 , and in the auxiliary problem 
F aux = F^ UX + F 2 UX , and T aux = T* UX + T 2 UX . Using these simplifications the general relation, 
Eq. (31), simplifies to 



(U aux • F a + ft aux • Ti) + Ui ■ F aux = jf (i + Bj ■ j^j : Vu 



dV- (32) 



We now need to find an auxiliary problem that facilitates the determination of the first 
order propulsion velocity, Ui, in the main problem. An appropriate candidate is the trans- 
lation of two touching spheres along the line of their centers without rotation, f2 aux = 0. 
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The exact analytical solution was given by Cooley and O'Neill 70 . By choosing this auxiliary 
problem, the relation further simplifies to 



-U a 



Fi + Ux ■ F a 



v f 



7o + B io ■ 7o : Vu 



dV- 



(33) 



If we do not allow the spheres to translate along the z-axis, Ui = 0, an external force, Fi, 
is required to hold the snowman in place given by 



■U aux • Fx 



Vf 



dV- 



(34) 



On the other hand, if we allow the snowman to translate freely without imposing any 
external forces, Fi = 0, then the first order propulsion velocity, Ui, can be determined from 



Ux F a 



Vf 



7o + B io ■ 7o : Vu 



dV, 



(35) 



where both F aux and the integral are expressed in terms of known Newtonian solutions of 
the main and auxiliary problems. Since the propulsion velocity Ui (with magnitude Ui) 
and the force in the auxiliary problem F aux (with magnitude -F aux ) act both vertically, the 
first order propulsion speed is finally given by 



aux JVf 



7 + Bj ■ 7 : 



dV, 



(36) 



where a positive value represents upward propulsion. 

Using Eq. (36) with the zeroth-order solution 61 and the auxiliary Newtonain solution 70 
we are able to determine theoretically the leading order propulsion speed of the snowman, 
U = De so Ui + 0(De so 2 ) = De(l — ()Ui + 0(De 2 ). The quadrature is performed in the 
tangent-sphere coordinates, with somewhat lengthy differential operations in evaluating the 
integrand. Our asymptotic results are shown in Fig. 2 as a solid line (blue online). We see 
that our results predict very well the propulsion speed of the snowman for small De when 
compared with numerical computations of the Oldroyd-B fluid (dot-dashed line - red online, 
in Fig. 2), and the agreement is excellent up to De ~ 1. 

Note that in order to compare the results between the second-order fluid calculation and 
the Oldroyd-B numerics, the dimensionless parameter B = — 2\E , 2 /^ / i i n the second-order 
fluid has to be taken to be zero because the second normals stress coefficient is zero in the 
Oldroyd-B model. Experimentally, indeed we have B < 1. Mathematically, the propulsion 
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velocity varies linearly with B, and a transition of propulsion direction occurs at B = 1. 
Such a transition also occurs in the direction of radial flow for a single rotating sphere in a 
second-order fluid (see Sec. V for a related discussion). 

C. Propulsion characteristics 

Anticipating the section where we make the link between rheology and locomotion, we 
now investigate the impact of the snowman geometry on its propulsion performance in 
the low-De regime where our asymptotic results via reciprocal theorem are quantitatively 
accurate. 

1. Touching spheres 

In the case of two touching sphere (h* = 1 + r*), the only free dimensionless geo- 
metric parameter is the ratio of the radius of the upper to that of the lower spheres 
r* = Rjj/Rl G [0, 1]. In the limit r* = 0, the snowman reduces to a single sphere, while the 
limit r* = 1 corresponds to two equal touching spheres; in both cases, there is no propulsion 
by symmetry. We therefore expect an optimal ratio r* for a maximum propulsion speed. 
Using the reciprocal theorem, Eq. (36), we calculate the propulsion speed as a function of 
r* (Fig 3, solid line - blue online) and compare with the numerical results in an Oldroyd-B 
fluid (Fig. 3, dots - red online) at De = 0.1 and a typical relative viscosity ( = 0.5. The 
asymptotic results agree very well with the Oldroyd-B computations. The optimal sphere 
size ratio occurs at r* pt « 0.58. In addition to our computations and theoretical calcula- 
tions, and based on physical understanding of the behavior of a single rotating sphere in a 
second-order fluid, a simplified analytical model can be constructed to predict the snowman 
dynamics with results shown as a dotted line (black online) in Fig. 3; the details of this 
simple model are given in Sec. V. 

2. Separated spheres 

Next, we let the two spheres be separated at a distance h* > 1 + r* (no longer touching). 
The two spheres still rotate at the same speed as a rigid body and the separation distance 
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FIG. 3. Propulsion speed of snowman with two touching spheres as a function of the ratio of radii 
at De = 0.1 and £ = 0.5. Dots (red online): numerical results in an Odroyd-B fluid. Solid line 
(blue online): theoretical calculation for second-order fluid. Dotted line (black online): simplified 
analytical model (Eq. 51). 

is kept fixed by connecting the spheres with a drag-less slender rigid rod (a mathemati- 
cally phantom rod) with negligible hydrodynamic contribution. Experimentally, this may 
be realized using, for example, using nanowires 71 . To compute the propulsion speed by the 
method described in Sec. IIIB2 and therefore Eq. (36), we need two new Newtonian solu- 
tions, namely the zeroth-order solution and the auxiliary problem. The zeroth order solution 
considers two separated unequal spheres rotating at the same rate in a Newtonian fluid, the 
exact solution of which was given by Jeffery 72 in bi-spherical coordinates. The appropriate 
auxiliary problem is the translation in a Newtonian fluid of the same two-sphere geometry 
along their axis of symmetry. Stimson and Jeffery 73 calculated that exact solution also in 
bi-spherical coordinates. 

For very separated spheres h* ^> 1, the propulsion is expected to decay with the sepa- 
ration distance. Hydrodynamic interactions between the two spheres is weak in this limit 
and each sphere behaves approximately as a single rotating sphere which does not propel. 
In Fig. 4a, the variation of the propulsion speed as a function of the separated distance 
is calculated for different fixed values of r*. The propulsion speed decays as expected for 
large h*. Interestingly, a non-monotonic variation occurs when the spheres are close to each 
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FIG. 4. Propulsion speed of a separated snowman as a function of (a) the separation distance, and 
(b) ratio of the radii, at De = 0.1 and £ = 0.5. Dots (red online): numerical results in an Odroyd-B 
fluid. Solid line (blue online): second-order fluid analytical calculation. Dotted line (black online): 
simplified model (Eq. 50). 



other (small h*). The swimming speed first increases with h*, reaching a maximum around 
h* ~ 2.5, before decaying to zero with further increase in h*. In Sec. V, a simple physical 
explanation to this non-monotonicity is discussed; the dotted line (black online) in Fig. 4a 
corresponds to the predictions by a simplified analytical model based on this explanation. 

For separated spheres, we can again vary the radii ratio, r*, at different fixed separated 
distance h* (Fig. 4b) and results similar to the case of touching spheres is observed: for any 
value of h* there exists an optimal value of r* at which the dimensionless propulsion speed 
reaches a maximum. The simplified model (Sec. V) again captures this trend qualitatively 
(dotted line - black online, Fig. 4b). 

Finally, by plotting the isovalues of the propulsion speed as a function of both r* and h* 
(Fig. 5), we are able to optimize the snowman geometry for the overall maximum propulsion 
speed. The optimal geometry occurs at (r*,h*) = (0.46,2.5), and a schematic diagram of 
the optimal snowman is drawn to-scale in Fig. 5. 
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FIG. 5. Optimization snowman propulsion. Iso-values of the dimensionless propulsion speed with 
dimensionless separation distance, h* , and ratio of sphere radii, r*. The optimal geometry for 
maximum propulsion speed is given by (r*, /i*) op t = (0.46, 2.5). A schematic diagram showing the 
optimal geometry is drawn to scale above. 

IV. MICRORHEOLOGY VIA SNOWMAN 
A. Scaling 

In the sections above we have derived an analytical expression, valid in the small De 
regime, relating the propulsion speed to the intrinsic properties of the complex fluid, namely 
the normal stress coefficients (Eq. 36). Turning all dimensionless variables back in dimen- 
sional form, this relationship reads formally 

U = (Cf*x + C?¥ 2 ) (37) 

where Cf and Cf are dimensionless coefficients depending solely on the snowman geometry 
(h* and r*) and defined as 

Cf = — ' (38) 

aux 

Q Svs (7S • To) : V*K n JV* 
CI = p- (39) 

aux 

Since the second normal stress coefficient ^2 is usually much smaller than the first normal 
stress coefficients ^1, we might ignore ^2 and obtain an estimation of \&i by measuring the 
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propulsion speed of a snowman U, i.e. 

*' " (40) 

where Cf depends only on geometry and can be computed using Eq. (38). This expression 
demonstrates the use of locomotion (U) to probe the local non-Newtonian properties of the 
fluid (tfi). 



B. Second Experiment: Repulsion of two equal spheres 



In scenarios where both values of \&i and \&2 are desired, a second experiment is neces- 
sary to obtain a second, independent, measurement of a combination of the normal stress 
coefficients. We propose to measure in the second experiment the relative speed (repulsion) 
of two rotating equal spheres of radius Re, with their centers separated by a distance h (see 
Fig. 6 inset for notations and geometry). Should the two equal spheres be connected as a 
rigid body, no propulsion would occur by symmetry. However, if the equal spheres are not 
connected but allowed to freely translate along their separation axis, upon imposing rotation 
they will translate with velocities of equal magnitude but opposite directions provided the 
fluid is non-Newtonian. 

We adopt the same non-dimensionalizations as previous sections (all lengths are now 
scaled by Re) and drop the stars for simplicity; all variables in this section are dimensionless 
unless otherwise stated. Denoting the dimensionless velocity of the lower sphere as V, 
we again expand the repulsion velocity in powers of De so , V = De^Vx + 0(De so 2 ), and 
determine the first order velocity Vi using our use of the reciprocal theorem as described in 
Sec. Ill B 2. By symmetry, the upper sphere translates with velocity —V (equal speed but 
opposite direction as the lower sphere). 

In this scenario we have to again define two setups, one for the main problem and one 
for the auxiliary problem. For the main problem, we consider the rotational motion of two 
free equal spheres about their line of centers 72 . Since the motion is force-free (F" = at 
each instant), Eq. (36) simplifies to 



2 

aux 



V* 



7 + B l0 ■ 70 : Vu au> 



dV, 



(41) 
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where fi" = for the same reason as explained in Sec. IIIB2. 

For the auxiliary problem, we consider the Newtonian translational motion (f2"ux = 0) 
of two equal spheres moving towards each other at the same speed and hence force, F^ ux = 
— F^ ux = F^ x . The exact solution to this problem was found by Brenner 74 in bi-spherical 
coordinates. We therefore have 



(Ul - U?) • Ff ux = jf (t o + Bj • 7 ) 



dV. (42) 



Note that the main problem here is a special case of that considered in Sec. IIIB2, but 
the auxiliary problem is completely different. We however still use the same symbols as in 
Sec. IIIB2 for simplicity. 

By symmetry, the two equal spheres propel with equal speed in opposite directions = 
— Uj = Vi , hence 

-2Vi • Ff ux = f 

Jv f L 

Since the repulsion velocity Vi (with magnitude V\) and the force in the auxiliary problem 
■^iux ( w hh magnitude i^ x ) both act vertically, the equation above can be rewritten as 



v 

7o + b 1q ■ 7o ) : Vu a 



dV. (43) 



Vi 



2F± 



V 

7o + B io ■ 7o ) : Vu a 



where a positive value of V\ represents repulsion. 



dV, (44) 



The only dimensionless parameter in this second experiment is the ratio of the separation 
distance to the radius of the spheres, which we write as h* = h/Ri- Using Eq. (44) we 
calculate the repulsion speed (V\ > 0) as a function of the dimensionless separation h* (solid 
line - blue online, Fig. 6, for De = 0.1 and ( = 0.5), and the results are found to be in 
excellent agreement with the Oldroyd-B calculations (dots - red online, Fig. 6). 

Back to dimensional variables, the leading order repulsion speed is formally given by 

V=(Cf* 1 + Ci* 2 )!^, (45) 

where Cf and Cf are dimensionless coefficients evaluated with the solution to the main and 
auxiliary problems described above 

c > = — IF • (46) 

aux 

w U* (To • To) : V*K nx dV* 
^ = ^ ^ (47) 

^^aux 
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FIG. 6. Dimensionless repulsion speed, V/ReQ, of two equal co-rotating spheres as a function of 
their dimensionless separation distance, h*, at De = 0.1 and £ = 0.5. Dots (red online): numerical 
results in an Odroyd-B fluid. Solid line (blue online): theoretical calculation for a second-order 



C. Determination of normal stress coefficients 

From measuring both the propulsion speed U of a snowman (given by Eq. 37) and re- 
pulsion speed V of the equal spheres (given by Eq. 45), we now have enough information to 
deduce both the first and second normal stress coefficeints (^i, ^2)- If we choose the same 
radius for the lower sphere in both experiments Re = Rl (we use Rl hereafter), we can 
write Eqs. (37) and (45) in a matrix form as 



where we denote by C the matrix containing the dimensionless coefficients (Cf , Cf , Cf , Cf ) 
in Eq. (48). The matrix C depends only on three geometric parameters, namely the ratio of 
the radii of the spheres in the snowman (r* = Ru/Rl), the dimensionless separation distance 
in the snowman (h^ = h/Ri) and that for the equal spheres in the second experiment 
(/ig = H/Re)- The coefficients of the matrix can be readily computed via Eqs. (38)-(39) 
and Eqs. (46)-(47), and thus the matrix in Eq. (48) can be inverted to obtain the values of 
v&i and \l>2- 



fluid. 




(48) 
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For practical implementation of this microrheological technique, measurement errors in 
the velocity of the snowman are inevitable and depend on the specific equipment employed 
for tracking the motion of the probe. However, the geometry of the snowman can be designed 
so that and ^2 are insensitive to measurement errors in the velocities U and V. The 
condition number (CN) of the matrix C to be inverted represents the maximum amplification 
factor of the relative measurement errors. The maximum relative errors in the normal stress 
coefficients would be equal to the condition number multiplied by the maximum relative 
measurement error. A small condition number is therefore desired. Similarly to the study 
by Khair and Squires 46 , we now investigate the value of condition number as a function 
of the geometry. To simplify the parametric studies, we first adopt the same separation 
distance in the first and second experiments (/ig = h% = h* = h/R^), and explore the 
dependence of the condition number on r* and h*, with results shown in Fig. 7a. The 
condition number does not vary monotonically with the parameters, which implies that 
optimization is possible. Under this requirement and within the ranges of values considered 
(r* G [0.2,0.98] and h* G [2.1,4]), the geometry yielding the lowest condition number is 
r* = 0.46 and h* = 2.1 (the corresponding condition number for C is ~ 27.6). When 
the requirement of h$ = is removed, by examining all combinations of the parametric 
values within the ranges (r* G [0.2,0.98], hg G [2.1,4], and h% G [2.1,4]), the minimum CN 
obtainable appears to be « 25.7 with r* = 0.46, h$ = 2.6, and h% = 2.1. 

The condition number can be further fine-tuned if we allow Re 7^ Re, in which case we 
have the new matrix relation 



The modified dimensionless matrix C in Eq. (49) now depends on one more parameter 
Rl/Re-i which is the ratio of the lower sphere radius in the snowman, R^, to that of the equal 
sphere, Re- In Fig. 7b, we investigate the dependence of the condition number of C with 
this new parameter, and adopt for all other parameters the optimal geometric parameters 
we determined before (r* = 0.46, hg = 2.6, and h% = 2.1). The variation turns out to be 
also non-monotonic, and a minimum is achieved when Rl/Re = 3.5 with CN 13.4. A 
schematic diagram showing the corresponding geometrical setup of the two sets of experiment 
is given to scale in the inset of Fig. 7. The condition number could be brought further down 
with a full four- dimensional parametric study and expanding the domains of the parametric 




(49) 
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FIG. 7. (a) Condition number (CN) for the matrix C as a function of sphere size ratio (r*) and 
dimensionless separation distance (h* = h$ = h^). (b) CN as a function of Rl/Re, for r* = 0.46, 
hg = 2.6, and h% = 2.1. 

studies. However, geometries yielding a lower CN may correspond to a negligible speeds 
undesirable for measurement. The current geometry (r* = 0.46, h$ = 2.6, = 2.1, and 
Rl/Re = 3.5) has both a relatively low CN and a high propulsion speed, making it ideal for 
experimental implementation. It is interesting to note that the optimal geometry for a small 
condition number we find here is close to the optimal geometry producing the maximum 
propulsion speed for the snowman (r* = 0.46, h$ = 2.5) determined in Sec. Ill C 2. 

V. QUALITATIVE PHYSICAL EXPLANATION 

In this section, we turn to an explanation of the physical origin of the non-Newtonian 
propulsion of a snowman. Based on physical intuition we present a simple model which 
successfully captures all the qualitative features of this mode of propulsion. 

We first look into the simplest related problem, that of a single sphere rotating in a 
complex fluid (a textbook problem discussed, for example in Ref. 9 ). Non- Newtonian stresses 
lead to the creation of a secondary flow in which the fluid moves towards the sphere in the 
equatorial plane and away from the sphere near the axis of rotation (see the inset of Fig. 8 
for an illustration of the secondary flow field) 9 . 

This secondary flow can be understood physically as a consequence of the hoop stresses 
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along the curved streamlines. Polymer molecules in the fluid get stretched by the flow, 
leading to an extra tension along streamlines. The presence of that extra tension along 
the closed circular streamlines leads to an inward radial contraction (like a stretched rubber 
band) pushing the fluid to thus go up vertically in both directions (by continuity) to produce 
the secondary flow. Notably, this secondary flow is independent of the direction of rotation 
of the sphere. 

The argument for locomotion of the snowman is then the following. Based on the one- 
sphere result, we see that when the two spheres in a snowman are aligned vertically and 
subject to a rotation they generate secondary flows and push against each other. For a 
single sphere, the strength of the secondary flow increases with the size of the sphere 9 . 
Consequently the smaller sphere is being pushed harder by the larger sphere than it is able 
to push against, and hence the two-sphere system is subject to a force imbalance, leading to 
propulsion. This physical understanding agrees with our results: propulsion always occurs 
in the direction of the smaller sphere, independently of the direction of rotation. Should the 
two spheres not be connected as a rigid body but free to translate vertically, they would 
repel each other, explaining physically our results in Sec. IV B. 

Based on this intuitive argument, we can now construct a simple mathematical model. 
Using the same notations as above, for a sphere of radius R rotating with an angular velocity 
O in a second-order fluid, the leading order solution v(<f),r,8) in spherical coordinates 9 is 
v* = v/RQ = (l/r*) 2 sin# + De so (l - £)[(l/2r* 2 - 3/2r* 4 + l/r* 5 )(3 cos 2 8 - 1) e r - 
3(l/r* 4 — 1/r* 5 ) sin# cos# eg] + O (De so 2 ). The Newtonian component of the flow field 
(De so = 0) is the primary flow field, and it has only a azimuthal ((f)) component. The 
secondary flow field, proportional to De so , is due to fluid elasticity and has only radial (r) 
and polar (8) components. As expected, the dimensional secondary flow v is quadratic in Q, 
confirming our physical intuition that it should be independent of the direction of rotation 
of the sphere. In the case where B < 1 (recall that B = —2^2/^1), the relevant limit for 
polymeric fluids, the secondary flow occurs in the direction intuited above and shown in the 
inset of Fig. 8. Note that the secondary flow field of a rotating single sphere would switch 
its direction when B went above one, explaining the switch in the propulsion direction of a 
snowman reported in Sec. IIIB2 in that limit. 

The dimensionless fluid velocity along the vertical axis (8 = 0) is given by v*(r, 8 = 0) = 
De so (1/r* 2 — 3/r* 4 + 2/r* 5 ) e r (8 = 0), where we have set B = to allow comparison with 
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FIG. 8. Rotation of a single sphere in a second order fluid. Radial velocity along the rotation axis 
(6 = 0) as a function of r/R at De = 0.1 and £ = 0.5. Inset: streamline pattern and velocity 
(shaded/color map) of the corresponding secondary flow. 

the numerical results. The velocity along the vertical axis, shown in Fig. 8, is expected 
to display non-monotonic variation with the distance from the sphere since the velocity 
decreases to zero both in the far field and on the solid surface. This is at the origin of the 
non-monotonic dependence of the snowman propulsion speed with the separation distance 
between the spheres shown in Fig. 4a. 

The forces experienced by the upper and lower spheres can be approximately estimated by 
considering the individual flow fields generated by their own rotation without the presence 
of the other sphere. We place an upper sphere at a distance h* = H/Rl from the center of 
the lower sphere. Using the same notations as in the previous sections, the dimensionless 
velocity generated by the lower sphere, and evaluated at the location of the upper sphere, 
is given by v* L (h*) = De so (1/h* 2 - 3/h* 4 + 2h* 5 ) e^(9 L = 0), where e^(9 L = 0) is the unit 
radial vector in the polar direction 9 L = 0, with respect to the coordinates system at the 
center of the lower sphere. Similarly, the dimensionless velocity generated by the upper 
sphere at the same distance, h*, but measured from the center of the upper sphere is given 
by v^(h*) = r*De so [{r*/h*f - 3 [r*/h*f + 2 [r*/h*f] e^{9 u = vr), where e^(9 u = tt) is 
the unit radial vector in the polar direction 9 U = ir, with respect to the coordinates system 
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at the center of the upper sphere. Note that e^{6 u = ir) = —e^(9 L = 0). As a simple 

approximation, we estimate the viscous drag force experienced by the upper and lower 

spheres to be ~ 67rr*v£(/i*) and ~ 67rv^(/i*) respectively. The difference between 

these two forces results in a net propulsive thrust. When dividing by an approximation of 

the translational resistance of the snowman at zero Deborah number with no hydrodynamic 

interactions, 67r(l + r*), we obtain a simple estimate of the dimensionless propulsion speed 

as U* « \F* V + F£ |/6tt(1 + r*). This leads to 

rr* |rX(^)+v^(^)l _ n r* [3fr*(r* 4 - 1) - h* 3 (r* 2 - 1) - 2r* 5 + 2] 

(1 + r*) s ° h* 5 (l + r*) ' 1 j 

In Eq. (50), we verify that U* vanishes when r* = (single sphere) and r* = 1 (equal 

spheres). For the case of touching spheres (h* — 1 + r*), Eq. (50) simplifies to 

U* ~De 2r * 3(1 ~ r * } (51) 

Does this simple model capture the essential propulsion characteristics? In Fig. 3, we plot 
the dimensionless propulsion speed of a touching snowman estimated by this simple model 
(Eq. 51) as a function of r* (dotted line - black online) and compare with the theoretical 
results from the reciprocal theorem approach (solid line - blue online) and the numerical com- 
putations (symbols - red online). The simple model correctly predicts the order of magnitude 
and captures qualitatively the variation with r*. For non-touching snowman, the qualita- 
tive model (Eq. 50) also captures qualitatively the variation of the dimensionless propulsion 
speed with h* (dotted line - black online - for r* = 0.6, Fig. 4a), also predicting an optimal 
separation distance and therefore supporting our understanding of a non-monotonic depen- 
dence with h* as arising from the non-monotonicity of the single-sphere velocity (Fig. 8). As 
expected, Eq. (50) also captures the non-monotonic variation with respect to r* for separated 
snowman (dotted line - black online - for h* = 8, Fig. 4b). 



VI. DISCUSSION AND CONCLUSIONS 

In this work, we present the design and mathematical modeling for a new non-Newtonian 
swimmer - the snowman - which propels only in complex fluids by exploiting asymmetry 
and the presence of normal stress differences under rotational actuation. The simple shape 
of our swimmer makes it ideally suited for experimental measurements. Note that if kept in 
place, the snowman would then act as a micro-pump for complex fluids. 
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The propulsion characteristics of the snowman are investigated by a combination of nu- 
merical computations (moderate values of De in an Odroyd-B fluid) and analytical treat- 
ment (small De in a second-order fluid). The underlying physics of propulsion, relying on 
elastic hoop stresses and geometrical asymmetry, is explained and based on this physical 
understanding a simple analytical model capturing all qualitative features is successfully 
constructed. Note that since, as a rule of thumb, inertial and elastic effects tend to produce 
secondary flows in opposite directions 9 , we expect that an inertial (instead of viscoelastic) 
snowman should swim in the opposite direction (from small to large sphere). 

The two-sphere setup proposed in this work is arguably the simplest geometry able to 
swim in a complex fluid under uniform rotation. It of course simplifies the analysis since 
the required Newtonian solutions to be used in our integral approach are all available. Any 
axisymmetric but top-down asymmetric geometry should also work, for example a cone, and 
clearly there remains room for shape optimization in that regard. Additionally, studying 
the snowman dynamics under a time-varying rotation could lead to a rich dynamics with 
potentially non-trivial stress relaxation effects. 

One of the main ideas put forward in this work is the use of locomotion as a proxy 
to probe the local non-Newtonian properties of the fluid. The snowman can be used as a 
micro-rheometer to estimate the first normal stress coefficient on its own, or to measure 
both the first and second normal stress coefficients with the help of another complementary 
experiment. Khair and Squires 46 recently proposed to measure normal stress coefficients 
by pulling microrheological probes and measuring the relative forces on the probes. In 
our work, we propose alternatively to perform only kinematic measurements of the sphere 
speeds instead of forces, which could present an interesting alternative from an experimental 
standpoint. 

We finally comment on a potential experimental implementation of the snowman tech- 
nique. We are aware of a number of rotational micro-manipulation techniques (see a short 
review in Ref. 75 ). For example, spinning micro-particles may be achieved by the use of 
optical tweezers and birefringent objects 76 . Birefringence allows the transfer of angular mo- 
mentum from the circularly polarized laser to the particle, producing controlled rotation. 
By rotating spherical birefringent crystals (vaterite), this technology has been implemented 
as a micro-viscometer to probe fluid viscosity 75 ' 77,78 . A similar mechanism may be useful 
for the two-sphere setup in this work although simultaneous rotation of two spheres may 
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introduce experimental challenges. Our dual-purpose snowman, both a micro-propeller and 
a micro-rheometer, invites experimental implementation and verification. 
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